Size effects in the exchange coupling between two electrons in quantum wire quantum 

dots 
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We theoretically investigate the properties of a two-electron system confined in the three- 
dimensional potential of coupled quantum dots formed in a quantum wire. For this purpose, we 
implement a variational Heitler-London method that minimize the system energies with respect 
to variational parameters in electron trial wavefunctions. We find that tunneling and exchange 
couplings exponentially decay with increasing inter-dot distance and inter-dot barrier height. In 
the quasi-one-dimensional limit achieved by reducing the wire diameter, we find that the overlap 
between the dots decreases, which results in a drop of the exchange coupling. We also discuss the 
validity of our variational Heitler-London method with respect to the model potential parameters, 
and compare our results with available experimental data to find good agreement between the two 
approaches. 

PACS numbers: 73.21. Hb, 73.21. La, 73.21.-b 



I. INTRODUCTION 

Very recently, the exploration of new hardware 
schemes at the frontier of solid state quantum computing 
has unveiled a new approach to fabricate coupled quan- 
tum dots (QDs) with controlling gate grid adjacent to an 
InAs quantum wire (QW)i In these device structures, 
electrons are laterally confined {i.e., perpendicular to the 
axial direction of the wire) by the wire external surfaces 
(wire diameters are tens or even a few nanometers),— and 
longitudinally confined in the wire axial direction by the 
electrostatic potential barriers created by the local con- 
trolling gates. The local gate width and separation range 
from ~ 10 to ~ 100 nm, which results in small effective 
dot sizes and inter-dot separations, so that size quantiza- 
tion effects and exchange coupling between the QDs are 
expected to be significantly larger than that in the two 
dimensional electron gas (2DEG) based semiconductor 
QDs. 3 In quantum wire quantum dot (QWQD) systems, 
the distance between the controlling gates and the QD 
region (~ 25 nm)i is smaller than that in 2DEG-based 
QDs (~ 100 nm)(i leading to better electrostatic con- 
trol of the charge (spin) states in the QDs. Furthermore, 
QWQD structures offer linear scalability (i.e., with the 
linear grid of the controlling gates) instead of the 2D 
scalability resulting from top or side gate patterning in 
2DEG-based QDs-i 

In laterally coupled 2DEG-based QDs, electron cou- 
pling occurs between the two QDs in the same plane 
as the 2DEG, and carrier confinement is much stronger 
in the perpendicular direction4££ In vertically coupled 
2DEG-based QDs, carrier confinement is weaker in the 
2DEG plane than that in the coupling (vertical) direc- 
tion (see e.g., Refs. 0and[||and references therein). The 
electron confinement and coupling defined in the fabrica- 



tion processes of coupled QWQDs considerably deviate 
from those achieved in the 2DEG-based coupled QDs: 
the electrons are strongly confined in the plane perpen- 
dicular to the axial direction of the wire because of the 
small wire diameter, while quantum mechanical coupling 
is achieved between two quantum wells with relatively 
weaker confinement, because of the controlling gate spac- 
ing and biases. 



While a wealth of literature has been dedicated 
to the theoretical study of 2DEG-based coupled 
QDs^iii^^^^^il less attention is paid to the 
QWQD systems. Among all investigated approaches, the 
Heitler-London (HL) technique is relatively simple in its 
conceptual methodology to extract the exchange coupling 
between coupled QDs&i£ its validity has been discussed 
for systems of various dimensions, 18 and efforts have been 
pursued to improve the energy calculation by integrating 
variational parameters in the HL methodi 19 ' 20 

In this paper, we compute the electronic structure of 
coupled QWQDs containing two electrons with a varia- 
tional Heitler-London (VHL) method. We first construct 
a three-dimensional (3D) model confinement potential for 
the QWQDs and introduce three variational parameters 
in the HL wavefunctions that account for the specific 3D 
confinement profile. We then numerically minimize the 
QWQD energies with respect to these parameters, and 
obtain the quantum mechanical and exchange couplings 
between the two electrons, as well as the addition energy 
of the second electron in the dot. In our analysis, special 
emphasis is placed on the geometric effects in the cou- 
pled QWQDs. We discuss the limitations of our VHL 
method but indicate its improvement over the conven- 
tional HL method. We finally compare our results with 
the available experimental data. 
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FIG. 1: (color online) (a) Schematic of coupled QDs Dl and 
D2 formed in a quantum wire. Gates Gl and G5 define the 
outer barriers of the two QDs; G3 controls the inter-dot cou- 
pling; G2 and G4 are plungers tuning the confinement in each 
QD. Charging current flows along the wire from source to 
drain, (b) Schematic of the confinement potential of the cou- 
pled QDs along the z (wire axial) direction. v?z,(r) and v 3 fl( r ) 
denote the localized s states in the left and right QDs, respec- 
tively. 



II. MODEL AND METHOD 

Figure [Ua) shows a schematic of coupled QDs Dl and 
D2 formed in a single quantum wire: gates Gl and G5 de- 
fine the outer barriers of the QDs, G3 controls the inter- 
dot coupling, and G2 and G4 are used as plunger gates 
for fine tuning of the potential in each QD. The charging 
current flows from source to drain along the wire. The 
material under consideration is InAs, for which we use the 
electron effective mass m = 0.023rno (Ref. |2~lI ) and dielec- 
tric constant e = 14.6. Hence, the effective Bohr radius 
ro = h 2 e/me 2 = 33.6 nm and effective Rydberg constant 
Ry = me 4 /2e 2 h 2 = 1.468 meV. We assume a parabolic 
confinement potential in the sy-plane V(p) = mil 2 p 2 /2, 
wherein we take ft p = H/m(D/2) 2 , and D is the nomi- 
nal value of the wire diameter. In the z-direction (along 
which the QDs are coupled), the confinement potential 
is modeled by a linear combination of three Gaussians: 



2m 



h i = Tr-[Ppi + ~Ai) +V( Pl ) + —p z j + V( Zl ), (4) 
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Note that we separate the motion of the electron in 
the xy-plane and in the z-direction in the single-particle 
Hamiltonian hi. In this work, we only consider mag- 
netic fields applied in the z-direction for which A = 
(—yBx + xBy) /2. Such a magnetic field effectively en- 
hances the confinement of the in-plane (xy-plane) ground 
state while preserving its cylindrical symmetry. 

In order to obtain the system energies, we use the fol- 
lowing trial wavefunctions: 



X±(r) 



(p L (r)±(p R (r) 



*± 0i,r 2 ) 



lL(ri)l,R(r2) ± ^L(r 2 )yfl(ri) 
V / 2(1±S 5 ) 
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(7) 



In above, X-h X-> an£ i ^- denote the single-particle 
ground and first excited states, two-electron singlet and 
triplet states, respectively. S —< ifiL\<pii > is the overlap 
between s orbitals ^i(r) and (pn(r) localized in the left 
and right QDs, respectively, and their specific expressions 
are 
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+Hcxp [ -p- ] , 

l bz 



(1) 



where Vq gives the depth of two Gaussian wells describing 
the confinement of the two individual QDs (we fix Vq = 
20 meV), Vb controls the barrier height between the two 
wells (Vb = except otherwise specified), l z is the radius 
of each QD, 2d is the nominal separation between the two 
QDs, and lbz denotes the radius of the central barrier. A 
schematic of V(z) is shown in Fig. HJb) by the solid line. 
The two electrons in the coupled QDs are described by 
the following Hamiltonian: 



Figure QJb) shows the schematic of <Pl( v ) and fR(r) 
in the z-direction by dashed lines on top of the poten- 
tial. With the variational wavefunctions, we calculate 
the single-particle ground and first excited state energies 



oO/l = 



energies 



X± 



E S/T 



X± ) i two-electron singlet and triplet state 



The detailed expres- 



sions of these matrix elements are given in the Appendix. 

In our VHL approach, we use the effective in-plane con- 
finement strength ui p , z-direction confinement strength 
uj z and effective half inter-dot separation a as variational 
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FIG. 2: (a) Single-particle ground state energy, (b) single- 
particle first excited state energy, (c) two-electron singlet state 
energy, (d) two-electron triplet state energy, (e) Coulomb en- 
ergy in the singlet state, (f) Coulomb energy in the triplet 
state, (g) tunneling coupling 2t, and (h) exchange interaction 
as a function of the half inter-dot separation d for l z = 30 nm 
and D = 20 nm. The inset in (g) shows 2t in the zoom-in 
region 20 < d < 30 nm. On each panel, the solid (dashed) 
line shows the VHL (HL) result. 



parameters to minimize the system energies^ By fix- 
ing these variational parameters equal to their nominal 

values lo = \l£l 2 p 



J% with lo c = eB/mc, lu z = Q z = 

y // 2Vo/ ml z and a = d, we recover the results from the 
conventional HL method. We calculate the Coulomb en- 
ergies in the singlet and triplet states by 



E, 



S/T 
Coul 



<*±|C|¥±> 

l±S 2 ((VHPR \C\ip L tp R ) 

± {<pL<fiR |C| <Pr<Pl)) , 



(9) 



(lPlPr \C\ VlVr) 

= (^ L (n)^(r 2 ) \C\ ^(n)^(r 2 )) , (10) 
(lPl^r \C\ <PrVl) 

= (<PL{ri)<PR{r 2 ) \C\ v3 fl (ri)^ L (r 2 )) . (11) 

The same notation has been used in expressing the ma- 
trix elements in the Appendix. Using both HL and VHL 
methods, we calculate the tunnel coupling 2t — e 1 — e° 
and the exchange coupling J = E T — E s . From the two 
electron wavefunctions, we compute the electron density 
as [ip L / R (r) are real] 



nS/T 



(n) 



2 y"|* ± (n,r 2 )| 2 dr 2 
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III. RESULTS 



(12) 



where C — e 2 /e|ri — r 2 |, and we have used the notation 



In Fig. [21 we plot (a) the single-particle ground state 
energy e°, (b) single-particle first excited state energy 
e 1 , (c) two-electron singlet state energy E s and (d) two- 
electron triplet state energy E T as a function of the half 
inter-dot separation d for l z = 30 nm and D = 20 nm. 
The solid and dashed lines show the results obtained 
from VHL and HL methods, respectively, from which we 
see that VHL method indeed gives lower system energies 
than the HL method. Here, we note that each energy is 
minimized with respect to a set of its own variational pa- 
rameters. We also note that the single-particle energies 
are positive simply because of the large energy contri- 
bution from the in-plane confinement: for D = 20 nm, 
hujp w 33 meV and is changed by less than 1% by varying 
d. 

For l z — 30 nm and d = 20 nm, the two Gaussian 
wells in Eq. fT]) are strongly coupled. As a result, the z- 
direction potential has a single minimum at z = 0, corre- 
sponding to a single QD. As d increases, a potential bar- 
rier between the QDs starts to emerge (for d > 21.2 nm). 
Meanwhile, the potential minimum is raised, and the 
z confinement in each individual QD becomes stronger. 
The behavior of the single-particle energies is a result of 
these combined effects. For example, as d increases from 
20 to 38 nm, both e° and e 1 sharply increase due to the 
large increase of the potential minimum [Figs. HJa) and 
(b)]. For 38 < d < 60 nm, e° still slowly increases, while 
e 1 starts to decrease. Our analysis based on the varia- 
tional parameters shows competing effects of the kinetic 
and potential energies in this region: for e , the kinetic 
energy increase dominates a slight drop of the potential 
energy, whereas for e 1 , the potential energy increase is 
offset by the drop in the kinetic energy. For very large d, 
both e° and e 1 approach a constant value (18.53 meV), 
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FIG. 3: (color online) Main panel: Exchange coupling J (solid 
curve) and tunnel coupling 2t (dashed curve) as a function of 
the effective barrier height V h . Values of half inter-dot sep- 
aration d corresponding to different V b e ^ values are shown on 
the upper horizontal scale. Inset: z-direction potential profile 
at V b & ^ values meV (red, solid), 5 meV (green, dashed) and 
10 meV (blue, dotted). Corresponding Vb values are —14.71 
meV, —1.16 meV and 6.65 meV, respectively. Values of other 
parameter are: D = 20 nm, d = l z = 30 nm, and l b — 30 nm. 




which corresponds to the limit of two decoupled quantum 
wells. 

The behavior of E s and E T [Figs. C^c) and (d)] re- 
sembles that of e° and e 1 , albeit a drop for d > 41 nm 
is observed for both quantities. The similarity implies 
that the single-particle energies are the dominant contri- 
butions to E s and E T , whereas the decrease of Coulomb 
energy with increasing d [Figs. [2](e) and (f)] has a minor 
influence. It is seen that at fixed d, the Coulomb inter- 
action is stronger in the singlet state, due to the larger 
overlap (S) in the two-electron wavefunction, which is a 
signature of the Pauli exclusion principle. 

In Figs. [H (g) and (h), we plot the tunnel coupling 2t 
and exchange coupling J as a function of d, respectively, 
both of which exhibit exponential decay with increasing 
d (strictly speaking, the decay is slightly slower than ex- 
ponential). In these figures, the solid (dashed) line corre- 
sponds to the VHL (HL) result. A much larger decrease 
of J (~ 10~ 8 ) than 2t (~ 10 -4 ) as d increases from 20 
to 60 nm agrees qualitatively with the Hubbard model 
J cx (2t) 2 /Ujj, assuming that the intra-dot Coulomb in- 
teraction Uh retains the same order of magnitude as d 
varies. Figs. [2] (g) and (h) show a large difference be- 
tween the tunnel and exchange couplings obtained by us- 
ing the HL and VHL methods, from which we notice that 
the HL method substantially underestimates the coupling 
between the two electrons,— especially for large inter-dot 
separations. For example, at d = 60 nm, the VHL result 
of 2t (J) is ~ 10 (~ 100) times of the HL result. 

The inset in Fig. [3] indicates that both the effective 



FIG. 4: (color online) (a) Exchange coupling J as a function 
of wire diameter D and half separation between the QDs d, 
which is set equal to QD radius l z (d — l z ). (b) J as a function 
of D for different d = l z values (shown in the figure). The 
J value on each curve is normalized to its value at D = 1 
nm. For (d = l x ) = 20, 30, 40, 50, 60 nm, J(D = lnm) = 
2.33 x 10~\ 2.47x 10~ 2 , 3.53 x 10" 3 , 1.37 x 10~ 3 , 4.81 x 10~ 4 
meV, respectively. 



barrier height V b e (i.e., the energy difference between 
the minima of the potential and its value at z = 0) and 
the distance between the two QDs (i.e., the distance be- 
tween the two minima of the potential) become larger 
as Vb is increased. Consequently, both 2t and J ex- 
hibit nearly exponential decaj*i£ with increasing V b ^ as 
shown in the main panel of Fig. [3J similar to the quasi- 
exponential drop of these two quantities with increasing 
QD separation 2d [cf. Figs. OJg) and (h)]. Again, we 
observe that J decays at a much faster rate than 2t. In 
experimental QWQD devices, the effective barrier height 
between the two QDs can be tuned by varying the central 
gate biasji and our analysis shows that the magnitude of 
the exchange coupling can be controlled by proper bias- 
ing the central gate as in 2DEG-based coupled QDs. 10 

Figure[3Ja) displays the exchange coupling J as a func- 
tion of both the wire diameter D and the half separation 
(d) between the two QDs. Here, we set d = l z noting 
that in experiments coupled QWQDs are defined on top 
of a linear gate grid with a particular periodicity, 1 which 
indicates that the effective QD size and inter-dot sepa- 
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FIG. 5: (color online) Electron density plot in the z-direction 
for (a) singlet and (b) triplet states at d = l z = 30 nm. In 
each figure, the density is plotted at D — 1 nm (red, solid), 
D — 10 nm (green, dashed-dotted) , D — 40 nm (blue, dashed) 
and D — 80 nm (black, dotted). For each D, the density is 
normalized to its peak value. 



ration are approximately the same. For the confinement 
potential given by Eq. fT]). this configuration leads to a 
constant effective barrier height of 5.68 meV, indepen- 
dent of the value of d = l z . The nominal confinement 
strength for a single Gaussian well (Vo = 20 meV) with 
l z = 15 nm (0.45r ) and 60 nm (1.78r ) is Ml z = 24.27 
and 6.07 meV, respectively. For a wire diameter D = 1 
nm, the nominal confinement is M7 P = 1.33 x 10 meV, 
which physically corresponds to the quasi- ID limit of the 
systems with aspect ratio (A p /A z = y/Q z /Q p ) < 0.05 for 
the investigated range of d — l z from 15 to 60 nm. In 
the opposite limit, where D = 80 nm, Kl p — 2.07 meV, 
the aspect ratio A p /A z > 1.71. At fixed D, J exhibits 
exponential decay with d = l z in Fig. Ufa), where it is 
also observed that J decreases with decreasing D at fixed 
d = l z . This trend is shown explicitly in Fig. QJb) for 
different d = l z . For comparison, the data on each curve 
are normalized to the value of J at D = 1 nm. At fixed 
d = l z , as D is decreased from 80 nm, J decreases, and 
the decreasing rate becomes larger as D approaches 1 
nm, which is the quasi — ID limit. The faster dropping 
rate of J near D = 1 nm is due to fl p cx l/D 2 , and the 
influence of the variation of tt p on J becomes stronger at 
smaller D (through the Coulomb interaction). Here, we 
note that although the general trend of J is to decrease as 
D is made smaller, the decreasing rates are much larger 
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FIG. 6: Variational parameters uj z (shown as huj z ), half sepa- 
ration a, the overlap S, and the Coulomb energies as a func- 
tion of D at d = l z — 30 nm. Left (right) panels are for the 
singlet (triplet) state. 



for intermediate d — l z values than for small or large 
values. 

These effects of the wire diameter variation on the ex- 
change coupling are rather unexpected as they show that 
J depends on the wire confinement perpendicular to the 
coupling direction. In fact, we find that the D variation 
not only changes u> p , but also induces significant changes 
in uj z and a, which minimize the singlet and triplet state 
energies. One can directly visualize such changes by in- 
specting the electron density variation with respect to 
the wire diameter. In Fig. [5j we plot the electron den- 
sity [Eq. [fT5]l] for different D values (d = l z = 30 nm) 
in (a) the singlet and (b) triplet states, respectively. For 
the singlet state, as D decreases, the separation between 
the two density peaks becomes larger, and the width of 
each peak becomes smaller. Consequently, the overlap 
between the two electrons is reduced. Similar effects are 
observed in the density of the triplet state to a less extent. 

In Fig. O we plot the D dependence of uj z and a on the 
top two rows. Both variational parameters increase as D 
is reduced, and the relative increase is more significant 
in the singlet state than the triplet state. As a conse- 
quence, the overlap S =< ^lI^-R > = ex P( — rnuj z a 2 /h) 
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FIG. 7: Main panel: Exchange coupling as a function of the 
magnetic field applied along the wire without the Zeeman 
effect for D — 20 nm, d = l z — 30 nm. Inset: same as main 
panel but with Zeeman effect. 



FIG. 8: (color online) Addition energy of the second electron 
E a (2) as a fucntion of wire diameter D and half inter-dot 
separation d = h. 



between the localized s states decreases with decreasing 
D in both states, and the relative decrease is larger in 
the singlet state [Fig. [BJ third row]. Despite this effect, 
the Coulomb interaction (Ecoui) becomes stronger with 
decreasing D [Fig. bottom row] for both states, which 
is due to the reduced size in the xy-plane. We also per- 
formed analysis for different d = l z and observed similar 
behavior as shown in Figs. [5] and [5] 

In general, the influence of the D variation on the ex- 
change energy results from the fact that the two electrons 
in the 3D QWQD system respond to the variation of a 
single external parameter by adapting all the variational 
parameters via the minimization of the system energy. 
The response varies depending upon the values of other 
fixed external parameters, which leads to the different 
decreasing rates observed in Fig. Hfb), for instance. 

The in-plane electron confinement can also be en- 
hanced by applying a magnetic field (B) along the wire 
without reducing the wire diameter. As with reducing D, 
J drops with increasing B as seen in Fig. [JJ main panel. 
The drop is nearly linear at large B, which is smaller 
than the drop rate when D approaches 1 nm [cf. Fig. 
SJb)]. This is because the in-plane effective (variational) 

confinement strength uj p « \J®% + w ? an d u c oc B, while 

il p oc 1/D 2 . It should be pointed out that the relatively 
small J drop in Fig. [7] is obtained in the absence of the 
Zeeman effect, and it is well known that unlike the small 
g factor in GaAs (g ~ —0.44), InAs QWQD has a much 
larger g factor (2 to 15.5)^ for which the Zeeman effect 
is dominant over the orbital effect in the J dependence 
on B. For example, the inset of Fig. [JJ shows that for 
g = 8,— the Zeeman effect totally smears out the orbital 
effect illustrated in the main panel of Fig. [Jj which leads 
to a negative J for B > 1.1 T. 

Because we model the confinement in the xy-plane by a 
two-dimensional harmonic oscillator potential, the single- 
particle levels in that plane are given by the Fock-Darwin 



spectrum, whereby the energy separation between the 
ground and first excited states decreases as B increases 
(in contrast, this separation increases with decreasing D). 
In our calculations in Fig. [Jj we take D — 20 nm and 
d — l z = 30 nm. At B = 10 T, the separation is 16.44 
meV, which is considerably larger than the sum of the 
single-particle energy separation in the ^-direction (2.11 
meV) and the Coulomb energy in the triplet state (1.91 
meV). This observation validates the assumptions of the 
HL method in which the wavefunctions are taken as linear 
combination of localized Gaussians separated in the z- 
direction, and only the ground state in the xy-plane is 
taken into account. 

Experimentally, the measurement of the addition en- 
ergy is frequently performed to probe the energy levels 
of the QD.— The addition energy of the N-th. electron is 
defined as E a (N) = fi(N) - (i(N - 1), where fi(N) is the 
chemical potential of an iV-electron QD. Within the VHL 
method, we are able to calculate the addition energy of 
the second electron as E a (2) = fi(2) - fi(l) = E s — 2e°, 
where E s and e° denote the singlet state energy and the 
single-particle ground state energy, respectively. We plot 
E a (2) as a function of the geometric parameters d = l z 
and D in Fig. [S] In general, as the QDs become larger in 
size (larger d = l z or D), the addition energy decreases, 
for both Coulomb interaction and size quantization ef- 
fects are reduced. We find (not shown) that at fixed 
d — l z and D, the Coulomb energy between the two elec- 
trons are uniformly smaller than E a (2), which is due to 
the size quantization effects in the coupled QWQDs. 

IV. DISCUSSIONS 

A. Limitation of the variational Heitler-London 
method 

As an inherent drawback of the HL method, our varia- 
tional scheme breaks down when the overlap between the 
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localized s states is large, which occurs for small inter- 
dot separations. For example, in our calculations of the 
system energies, the VHL method fails for (d = l z ) < 12 
nm (0.57ro) independent of D. A signature of the VHL 
approach breakdown at small d is that the variational 
parameter a becomes zero in the minimization process. 
This numerical behavior stems from the fact that at small 
d a global minimum in the system energies does not exist 
for the physical range of a, given the expression of the 
variational wavefunction. We note that this shortcoming 
in the HL method in not apparent in the conventional HL 
approach. As long as (d = l z ) > 0, one can still use the 
HL method (without variation) to calculate the system 
energies even though the obtained result is likely to be 
unphysical. 

In Ref. Il8|, it was pointed out that the HL method 
breaks down as the quantity c = \J it / 2(e 2 / eas) / tv-Oo 
(a,B = ^/mwo) is larger than 1.95, 2.8, and 5.8 for 
coupled QDs with harmonic oscillator confinement Hluq 
in each direction for ID, 2D and 3D potential models, 
respectively (this is an extension of the result in Ref. 0). 
We investigate l z from 15 to 60 nm, which corresponds 
to c ranging from 0.44 to 0.87, and is uniformly smaller 
than the smallest breakdown value c = 1.95. However, 
as a check of this criterion, we extend our calculation 
to very large value of d = l z and find that for D = 20 
nm, J becomes very noisy and oscillates randomly for 
(d = l z ) > 206 nm, for which the variational parameter 
Huj z is 1.553 meV, corresponding to c = 1.723, which is 
similar to the ID limit claimed above. However, at this 
point, J ~ 10~ 14 meV, which bears no practical interest. 



QWQD system than those calculated by conventional HL 
method with the relative difference in the tunnel and ex- 
change couplings exceeding 100%. As in coupled GaAs 
QDs based on 2DEG, tunnel and exchange couplings ex- 
hibit exponential decay with increasing inter-dot distance 
or barrier height. Due to the 3D nature of the system, 
increasing the confinement in the in-plane directions re- 
duces the overlap of the two electrons in the coupling 
direction (along the wire), which results in the decrease 
of the exchange coupling. For QDs with different sizes, 
the addition energy of the second electron is found to 
be uniformly larger than the two-electron Coulomb in- 
teraction because of size quantization effects. By fitting 
the model potential to experimental parameters, we ob- 
tain exchange coupling in agreement with experimental 
data. Experimental structures based on InAs QWQDs 
may benefit from the relatively large exchange coupling 
towards quantum computing applications. 
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B. Comparison with experiments 

In recent experiments on InAs QWQDs, J = 2.8 
to 3.2 meV was reported for a single QD formed in 
a wire with effective harmonic confinement strength 
Ml z = 6.3 meV (corresponding to confinement length 
2\ z = 2^h/mtt z = 46 nm) and Ml p = 40 meV 
(2A P = 2y/h/m£l p — 18 nm).— By fitting these values 
in our model (D = 18 nm, Vq = 41.6 meV, Vb = meV, 
d = nm and l z = 117.9 nm), we obtain J = 3.51 meV, 
which is comparable to the experimental result. 

We note that J ~ 3 meV as obtained above is the result 
for a single QD with potential minimum at z — 0.— For 
double QDs with D = 20 and d = l z = 30 nm, we obtain 
J ~ 0.5 meV (Fig. [7]), which corresponds to a time scale 
(tj = h/J) ~ 1.3 ps, on the same order as the reported 
spin decoherence time T 2 = 0.5 — 1 ps in InAs QWQDs 21 
and much smaller than the reported spin dephasing time 
T 2 * = 50 - 500 ps in self-assembled InAs QDs^ 

V. CONCLUSION 

By introducing variational parameters in the HL trial 
wavef unctions, we achieved lower energies of coupled 



APPENDIX 

The single particle Hamiltonian h can be rewritten as 



h = h° L/R + W L/R , (A.l) 
*°/* = 2^( Pp + 1 c A ) 2+1 2 m "> 2 

W L/R = \m (n 2 p - w%) p 2 + V(z) - ~mwl(z ± a) 2 . 

(A.3) 

Since h° L/R ip L/R (r) = E tp L/R (r) and E = huj p + huj z /2, 
we only need to calculate the matrix element of Wl/ r . 
Thus, we have 
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e 0/1 = (X± 



X± 



1±S 



+v b 



mu) z l 2 



+ 1 



exp 



1 raw. 



1 exp 



1 raw 



itilvJ 2 



{a-df 



exp 



1 raw. 



mivj 2 



il n 



moj z 



(a + df 



exp 



1 raw. 



11 h 



mi 



± 5H 



S 



2„2 



FT" + 1 T — muj z a 

muj z V hz 



(A.4) 



In above, S = J dr(p* L (r)ipn(r) = exp (— mLo z a 2 /h) is the The singlet and triplet energies are evaluated in a sim- 
overlap between the two localized s states. ilar fashion and the results are: 



E S/T 



= 2E + 



H orb 

1 



1±S 2 



[(<Pl<Pr\Wi +W 2 + W 3 + C\ip L <p R ) ± Re (<p L <p R \Wi + W 2 + W 3 + C\ip R <p L )] , (A.5) 



(VL^R \V{ z l) + V(z 2 )\ VlVr) 

h 



-2V 
+ exp 
+2V b 



mu) z l 2 



1 



exp 



I mu> 2 
1 

II h 



1 muj z x mu z 5 



mw z l 2 z 



exp 



1 raw. 



(</?l</?r I Wi| ^flV'i} = (^i^ii \V(zi) + V(z2)\ VrVl) 



= 2exp^- 



-2V 



m,Lu z l 2 



+ 1 



muj z 



(a + d) 2 



«2, 



z 2 

-a 



exp 



1 mw ; 



mi 



(A.6) 



+ v b 



mw z l 2 z 



(A.7) 
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(ifiL^R \W 2 \ ifL<Pll) 



{VLfR \W 2 \tp R ip L ) 



(ifiL^PR \W 3 \ip L (f R } 



VlVr 
-mu; 2 z 

h 



-~mw 2 z {zi + a) 2 



-muj 2 (z 2 - af 



1 



--mwl{z x + a) 2 - ^mLu 2 (z 2 - af 



2muj 
1 



+ a exp 



Jz „2 



mini 



to 



po) [Pi + P2) 



2 + ^ 



(ipL<PR |W 3 | VrVl) = (fLfR 



-m (n 2 p - u> 2 p0 ) {p\ + p 2 2 ) 



VlVr 



fR^L 



— (n 2 + 



j^J exp (- 



VlVr 



<Pr<Pl 



(A. 



(A.9) 



(A.10) 



(All) 



The Coulomb matrix elements in Eq. ([9]) are given by 



(<Pl<Pr \C\ (Pl<Pr) = ( PLVR 



(<Pl<Pr \C\ <Pr<Pl) 



e ri - r 2 



VlVr 



e 2 / 2mu>, 



VlVr 



o 1 + 



— exp ( -2— -a 2 t 2 ) dt, 
f2 V n 



e |ri - r 2 

e 2 / 2mw 



e \ nti J 6XP (" 2 — ° ) 



e 2 (2muj z \ 2 



o 1 



It 2 



exp j x < 



1 _ a* 



arctanh 



1 arctan 



1 _ <^ 



(A.12) 



for w 2 < 



for oj 2 > lu 



for w 2 = lu„ 



(A.13) 



r 



Practically, the one dimensional integrals in Eqs. (A.12) 
and (|A.13|) ] are numerically evaluated using adaptive 
quadratures. We note that in the ID limit (oj p — > 
oo), the integrals have logarithmic divergence j*S while 
they both approach zero in the oppsite limit (lo p — » 



0). 



For w„ 



ojq, the integrals in Eqs. 



(A.12) and KJg| simpl ify to [e 2 /(2ae)]Erf(a v /2mtj /?i) 
and {e 2 /e)y/ 2mwo / (7rS) ) exp ( — 2mcuoa 2 /H) , respectively. 
These results are identical to the results in Ref. LL8|, where 
the Coulomb matrix elements were calculated between 
coupled spherically symmetric Gaussian trial wavefunc- 
tions. 
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